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ABSTRACT 

The eclipsing low-mass X-ray binary 4U 1822-371 is the prototypical accre- 
tion disk corona (ADC) system. We have obtained new time-resolved UV spec- 
troscopy of 4U 1822-371 with the Advanced Camera for Surveys/Solar Blind 
Channel on the Hubble Space Telescope and new V- and J-band photometry 
with the 1.3-m SMARTS telescope at Cerro Tololo Inter- American Observatory. 
We use the new data to construct its UV/optical spectral energy distribution 
and its orbital light curve in the UV, V, and J bands. We derive an improved 
ephemeris for the optical eclipses and confirm that the orbital period is changing 
rapidly, indicating extremely high rates of mass flow in the system; and we show 
that the accretion disk in the system has a strong wind with projected velocities 
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up to 4000 km s _1 . We show that the disk has a vertically-extended, optically- 
thick component at optical wavelengths. This component extends almost to the 
edge of the disk and has a height equal to ~0.5 of the disk radius. As it has a low 
brightness temperature, we identify it as the optically-thick base of a disk wind, 
not as the optical counterpart of the ADC. Like previous models of 4U 1822-371, 
ours needs a tall obscuring wall near the edge of the accretion disk, but we in- 
terpret the wall as a layer of cooler material at the base of the disk wind, not as 
a tall, luminous disk rim. 

Subject headings: binaries: eclipsing - pulsars: individual (4U 1822-371) 



Introduction 



4U 1 822-371 is a low- m ass X-ray binary star (LXMB) with an orbital period of P or b = 
5.57 hr (IMason et ajj Il980l ). The compact star in the sy stem is an X-ray pulsar with 
spin period P pu i S e = 0.593 s (IJonker fc van der Klia 120011 ) and is, therefore, a neutron 
star. The UV7 optical spectrum of 4U 1822-371 is rich with emission and absorption lines 



1982 



2003: 



Jonker et al. 



(Charle s et al.lll 980: Maso n fe Cordovall 198 2a: Ma son et al.lll982tlHarlaftis et al.lll997l : ICowlev et al 



20031 : iHutchings et al. 



20051 ). Doppler tomography of the He II A4686 



and O VI A3811 emission lines confirms expectations that there i s an accretion disk around 
the neutron star fed by a stream of gas from the secondary star ( jCasares et al.l 120031 ). 



4U 1822-371 is one of the rare eclipsing LXMBs (IMason et a 



Many analyses of the ecl ip se have been published (IWhite et al. 



Mason fc Cordovalll982al Jbl; iHellier fc Masonlll989l : iBaptista et al. 



1980: 



981 



White et al 



White fc Holt 



2002 



Cowley et al. 



1981). 



1982 



2003|) 



and, although the detailed results of the various analyses often differ greatly, there is sub- 
stantial agreement on several points. The eclipse is a transit of the secondary star across the 
accretion disk and, since the eclipse is very bro ad at optical wavelengths, t he accretion disk is 
large, altho ugh exactly how large is uncertain: IMason fc Cordova! (Il982bl ) find r disk /a = 0.28 
to 0.38 and IHellier &: Masonl (119891 ) find r^ s ^/a = 0.58 ± 0.08, where r^sk is the radius of 
the disk and a is the separation of the stars. The X-ray eclipse is also broad. Its flux at 
minimum is at about 50% of the uneclipsed flux level and it has no sharp features, showing 
that the X-ray emission comes from a p artially-eclipsed ex t ended cloud around the n eutron 
star - an accretion disk corona (ADC) (IWhite et al.lll98ll ; IMason fc Cordova! Il982al ). The 
corona extends to ~ 0.5r disk in the plane of the disk and is vertically extende d. Estimates of 



the r a tio of its height to its radius in the plane of the di sk range from 0.3 to 1 (IWhite &: Holt 



19821 : IMason fc Cordova! Il982al : IHellier fc Masonl Il989f ). The orbital inclination lies in the 

range 78° < i < 84°. 
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The ratio of the X-ray luminosity of 4U1822-371 to its optical luminosity, Lx/L opt , 
lies between ~ 15 and ~ 65, much less than that for typical LXMBs, indicating that the 
ADC is optic ally thick and blocks most of the X-ra y emission from the inner disk and 
neutron star ([Griffiths et al.l Il978l ; iMason et al.lll980l ). Fits to the X-ray spectrum using 
realistic Comptonization models yield electron temperatures near T e pa 5 x 10 7 K and electron 
scattering optical depths in th e range r R = 13-26 depending on th e geometry of the ADC 
and the Comptonization model (jParmar et al.ll2000l ; llaria et al. The mean unabsorbed 



apparent luminosity of 4U 1822-371 in the 0.1-100 keV energy range is 1.15 x 10 erg s 



for a distance of 2.5 kpc (llaria et al.l l200ll ). As most of the X-ray flux escapes roughly 
perpendicularly to the orbital plane and is beamed away from the Earth, the true luminosity 
is much greater than the apparent luminosity. If one corrects for the ADC obscuration by 
simply multiplying L x by a factor to raise L x /L opt to the typical value for LXMBs (~ 500), 
the true X-ray luminosity is ~ 10 37 erg s _1 . Even this should be taken as a lower limit 
because the disk is observed edge-on, reducing L opt . A true luminosity near the Eddington 
luminosity is not impossible. 

All solutions for the X-ray eclipse have invoked a vertically-extend ed wall of optically - 
thick material whose height varies with angle around the neutron star (IWhite et al.lll9811). 



he w all generally has been identified as a vertically extended disk rim, and, i ndeed, 



Hellier fe Mason 



1989) showed the wall must be near the edge of the disk. The rim invoked by lHellier Sz Mason 
19891) has 0.08 < h r iml r &sk < 0.22, while the m aximum he i ght of the rim invoked by 
White fc Holtl (119821 ) approached h rim /r disk pa 0.5. llaria et al.l (120011 ) also needed a thick 



belt of absorbing material around the ADC to match the observed X-ray spectral energy 
distribution (SED). Placing the belt at the edge of the disk, they found the angle subtended 
by the absorbing region is 16°, or h r i m /rsak ~ 0.29. These rims are an order of magnitude 
taller than th e expe cted thickness of the disk proper and cannot be pressure supported. 
White fc Holtl (J1982J) invoked turbulence to support the rim but the turbulence would have 
to be highly supersonic, needing Mach numbers between 7 and 15. While it is not inconceiv- 
able that supersonic turbulence could be generated by the impact of the accretion stream 
on the disk, it is unclear how supersonic turbulence could be maintained around the entire 
circumference of the disk. 

It has been universally assumed that the only contributors to the UV/optical continuum 
are 1) the accretion disk, 2) its extended rim, which may be irradiated and have a higher 
temperature on its surface towards the neutron s t ar, and 3) the irradiated s econdary star 
(IMason fc Cordova! Il982bl ; iHellier fc Masonl Il989l ; iPuchnarewicz et al.l Il995l ) . It also has 
been assumed that the ADC is optically thin at UV/optical wavelengths, even though it is 
optically thick to electron scattering at X-ray wavelengths. These contradictory preclude a 
UV/optical contribution from the ADC or from vertical structures other than the rim. 



-4- 



The ADC and the eclipses lend particular importance to 4U 1822-371. ADC sources 
are more numerous than suggested by the small number that have been identified. The 
luminosity of an ADC is typically ~ 100 times less than the luminosity of the inner disk 
and neutron star, so ADCs are visible only in those systems with orbital inclinations near 
enough to 90° that the outer accretion disk or t he ADC itself block s most of the flux from 



the more luminous components of the binary (IWhite fc Holtl Il982l ). Many ADC sources 
must be lurking among the low-inclination LXMBs, perhaps 5-10 times more than are 
currently known. 4U 1822-371 is the prototype for these systems. Models for accretion 
onto neutron stars and black holes in X-ray binaries or onto supermassive black holes in 
AGN need much more than flat, unadorned accretion disks. They must also include high- 



(Naravan et al. 


1998 


2000 




Lasota 


2001 


). 1 



(INarayan et al.l 11998c iBlandford fc Begelmanl Il999l ; IWijers &: Pringld Il999i ; iNarayan et al. 



20011 ). ADCs are examples of high-temperature accretion flows and vertically- 
extended structures. The eclipses of 4U 1822-371 offer the opportunity to map these struc- 
tures in detail. 

We have obtained new EST objective-prism ultraviolet spectroscopy and V and J pho- 
tometry of 4U 1822-371, from which we construct its SED and its orbital light curves in 
the UV, V, and J bands (Section 2). In Section 3 we derive an improved ephemeris for the 
optical eclipses and confirm that the orbital period is changing rapidly, indicating extremely 
high rates of mass flow in the system. In Section 4 we show that the accretion disk has a 
strong wind with projected velocities up to 4000 km s _1 . The core of this paper is Section 5, 
in which we model the orbital light curves with our light curve synthesis program. We show 
that much of the disk is vertically extended and optically thick at UV/optical wavelengths. 
This vertical structure extends almost to the edge of the disk and has a height equal to ~0.5 
of the disk radius. As this structure has a brightness temperature near 3 x 10 4 K, we identify 
it as the optically-thick base of the disk wind, not as the optical counterpart of the ADC. 
We, too, need a tall obscuring wall near the edge of the accretion disk, but we interpret the 
wall as layer of cooler material at the base of the disk wind, not as a tall disk rim. Our 
results are summarized in Section 6. 



Observations and Data Reduction 



We observed 4U 1822-371 with the ANDICAM optical/IR phot ometer on the SMA RTS 
1.3-m telescope at CTIO using Johnson V and CTIO/CIT J filters JDePov et al.ll2003h . We 
obtained one or two observations per night over two observing seasons, the first season from 
2005 June 30 to September 15, the second from 2006 March 29 to October 20. The data 
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were reduced to instrumental magnitudes with IRAF0, and the instrumental magnitudes 
were converted to absolute fluxes using the photometric standard stars PG1657+078 for the 
V band and 2MASS J18254777-3706131 for the J band JLandoltlll992l ). both of which were 
observed on three photometric nights in 2007 June. 

We also obtained a time series of objective-prism ultraviolet spectrograms of 4U 1822- 
371 with the Hubble Space Telescope (HST) using the Advanced Camera fo r Surveys/Solar 
Blind Channel (ACS/SBC) and the PR130L prism JPavlovskv et al.l l2006h . We obtained 
300 spectrograms in two visits between HJD 2453828.7 and HJD 2453830.9 (2006 April), 
each visit lasting five HST orbits. The spectrograms were integrated for 45 seconds and 
separated by a 40-second dead time, with occasional long exposures of 350 seconds when 
data were transfered from the ACS internal buffer memory to the HST solid state data 
recorder. The spectrograms cover the wavelength range from 1222 A to 2002 A, but we 
used only the wavelengths from 1222 to 1900 A because there is a significant red leak in the 
SBC at wavelengths greater than 1900 A. The spectral resolution is R = A/AA ~ 150 near 
1300 A, but decreases to R « 100 near the C IV doublet and to R ps 40 near 1900 A. The 



UV data were reduced with the aXe package in STSDAS (IKuemmel et al.l 12004 ). 



The time-averaged mean UV spectrum of 4U 1822-371 is shown in Figure [TJ The 
prominent emission lines are N V at 1240 A, a blend of O IV, O V, and Si IV near 1370 A, 
the C IV doublet at 1548/1550 A, and He II at 1640 A. There are also interstellar absorption 
lines from blends of Si II/C I near 1260 A and O I/Si III near 1300 A. 



3. The UV/Optical Ephemeris and Orbital Light Curves 



Because the V and J light curves are sparsely sampled, it was not possible to extract 
times of minima for individual eclipses. Instead we measured a single m ean time of minimum 
for each season by folding the V^-band data on the orbital period from iParmar et al.l (120001 ) 
to form seasonal mean orbital light curves. The times of minimum light extracted from 
the mean light curves are T min = HJD 2453618.677 ± 0.005 for the 2005 data and T min = 
HJD 2453932.720 ± 0.004 for the 2006 data. To form UV light curves the flux in each 
UV spectrogram was integrated fr om 1222 A to 1900 A and then divided by AA = 678 A. 
We folded the light curves on the IParmar et al.l (120001 ) period and extracted two times of 



minimum light, one for each HST visit: T n 



HJD 2453828.972 ± 0.001 and T„ 



1 IRAF is distributed by the National Optical Astronomy Observatory which is operated by the Associ- 
ation of Universities for Research in Astronomy Inc., under the cooperative agreement with the National 
Science Foundation 
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HJD 2453830.830 ±0.001. 



The four new eclipse times plus all of the previously published optical eclipse times are 
listed in Table [TJ A fit to the eclipse times yielded the ephemeris 

T min = HJD 2445615.31166(74) + 0.232108641 (80) E + 2.46(21) x 10~ n £ 2 , (1) 

where E is the eclipse number. The O — C diagram for this ephemeris is shown in Figure [2J 
The quadratic te rm in the optica l ephem eris is consistent with the low-significance quadratic 
term derived by iBaptista et al.l (120021 ) and agrees with the qua dratic term in the X-ray 
ephemeris, (2.06 ± 23) x 10™ 11 , to within the measurement errors (IParmar et al.ll2000f) . The 



T m i n o f our new optical ephemeris lags the T min from the X-ray ephemeris of IParmar et al. 



(120001 ) by 100 ± 65 seconds. Although the lag is only marginally statistically significant, it 
cannot be dismissed because the X-ray and optical/UV flux emanate from different regions 
and could be eclipsed at slightly different times. Recently, Burderi et al. (private communi- 
cation) derived a new X-ray ephemeris for 4U 1822-371. Using this new ephemeris the lag 
is 122 seconds. 

The rate of change of the orbital period is P = 2.12 x 10 -10 and the timescale for a 
change in orbital period is P/P = (3.0 ±0.3) x 10 6 yr. One must generally be cautious about 
interpreting timescales for orbital period changes as evolutionary timescales because mass- 
transfer binaries often show short-term systematic advances and delays in times of eclipse 
that can mimic ra pid period changes. The X-ray binary EXO 0748-676 shows this behavior 



flWo 



ff et al.l 120091 ) and there are many examples among cataclysmic variables (see IWarner 



19951 ). However, the eclipse times for 4U 1822-371 span 27 years and the accumulated phase 
shift is A0 = AT m ; n /P « 0.17, an order of magnitude larger than the random phase shifts 
observed in EXO 0748-676 and cataclysmic variables. It appears safe, therefore, to identify 
the value of P/P measured from the ephemeris as the time scale for orbital evolution. 

On dimensional grounds one has m/m ~ \P/P\. If mas s and angular momentum are 
conserved, for example, then rriNs/i^NS = [3(1 — q)/q]P/P (IWarnerlll995l ). where mjvs is 
the mass of the neutron star, q = m^jvn^s^ an d Tn 2 is the mass of the secondary star. For 
q = 0.3 and m^s — 1.35M Q , we find m^s = 6.4 x 10~ 8 M Q yr -1 . This is an extremely high 
rate of accretion, co rresponding to <y 4 tim es the Eddington luminosity. In this we agree 
with the analysis of iHeinz fc Nowata (120011 ). We do not agree that the flow of mass and 
angular momentum can be characterized by a single parameter. Specifically, the angular 
momentum per unit mass can vary enormously in various parts of the flow and depends on 
the detailed properties of the flow. Absent a precise description for the flow of mass and 
angular momentum, it is not possible to infer accurate rates of mass transfer and mass loss 
from the evolutionary timescale, nor even a ratio of the tw o rates. Few c ataclysmic variables 
have rates of mass transfer greater than 3 x 10~ 8 M & yr -1 (jWarnerlll995l ). so we suspect that 
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the rate of mass transfer in 4U 1822-371 is at least a factor of two smaller than the value for 
conservative mass transfer. Yet, a mass transfer rate of 3 x 10 -8 M & yr -1 is still extremely 
high and much of the transferred mass cannot be accreting onto the neutron star, but leaves 
the system altogether. 

Figure [3] shows the UV light curves of 4U 1822-371 from the two HST visits folded on 
the orbital period given by equation [TJ and Figures 0] and [5] show the two seasons of V- 
and J-band photometry also folded on the orbital period. The folded UV light curves and 
the l^-band light curve from 2006 have been averaged into 200 equal-width phase bins. The 
J-band light curves from both 2005 and 2006 and the l^-band light curve from 2005 have 
fewer observations, and have not been binned. We define the internal error of a data point 
to be the standard deviation that includes photon counting noise and calibration errors. The 
internal errors of the UV fluxes are small, typically only 0.5%; the internal errors of V- and 
J-band data were typically 4 - 5% but occasionally somewhat larger. The error bars in the 
figures correspond to these internal errors. 

The light curve of 4U 1822-371 shows large-amplitude non-repeating variability ("flick- 
ering" ) and, because of this variabili ty, individual orbital light curves can differ substantially 



from the mean orbital light curve (IBaptista et al.l 120021 ) . Let us define the external error 



of a data point as the standard deviation introduced by the non-repeating variability. We 
estimated the external error of the UV fluxes by comparing the light curves of 4U 1822-371 
from the two HST visits. After folding all the UV data together to form a single (noisy) light 
curve, we divided the light curve into four phase sections (0.11 - 0.23, 0.33 - 0.52, 0.62 - 
0.9, and 0.9 - 1.11) and fit a polynomial to each section. The standard deviation of the data 
points about the fitted polynomial ranged from 5% to 9%. We equate the external error of 
the data points in each section to the standard deviation for the section. Much the same 
method was used to determine the external error in the V- and J-band data except that the 
light curves were divided into just two sections covering phases 0.25 - 0.70 and phases 0.70 - 
1.25. The error bars in Figure [TT1 correspond to these external errors. When calculating the 
X 2 of models fitted to the light curves (Section 5), we will use the larger of the internal and 
external errors for each data point. With the exception of a few data points in the V- and 
J-bands, the external errors are much larger than the internal errors and it is these large 
external errors that dominate the fit. 

The shape of the V- and J-band light curves changed systematically between 2005 and 
2006. While the 2006 light curves are roughly symmetrical about phase 0.0, the 2005 light 
curves are asymmetric, having a large hump near phas e 0.2. The asymme t ric light curves 



from 2 005 a re similar to the optica l light curves shown in lMason et al.l (119801 ) jHarlaftis et al. 



(119971 ). and IBaptista et al.l (120021 ). The more symmetric light curves from 2006 are similar 
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to our UV light cur ves - also obtained in 2006 - and to the optical light curves shown in 



Cowlev et al.l (120031 ). 



Figure [6] shows the X-ray light curve of 4U 1822-371 between 1996 January 5 and 2008 
March 5 from the All Sky Monitor (ASM) on the Rossi X-ray Timing Explorer (RXTEjE 
Each point in the figure is the monthly average count rate summed over the three ASM 
passbands. The dates on which optical eclipses were observed are marked on the figure. 
The X-ray light curve shows a low-amplitude variability superimposed on a slow decrease 
in flux. The pronounced changes in the optical/UV light curves from year to year were not 
accompanied by large changes in the X-ray light curve. 



4. The UV Spectrum, UV/Optical Spectral Energy Distributions, and 
Evidence for a Strong High- Velocity Wind 

Although the direction to 4U 1822-371 is within 11° of the direct ion to the galactic 
center, 4U 1822-371 is close to Baade's window and its reddening is low. iMason fc Cordova 
()1982al ) find a color excess E(B — V) = 0.1 with 3<r upper and lower bounds of 0. 29 and 0.01 



based on the depth of the 2200 A interstellar absorption feature. According to llaria et al. 



(120011 ) the photoelectric absorption at soft X-ray energies corresponds to a neutral hydrogen 
column density of N H g (1.05 ± 0.4 5) x 10 21 cm" 2 . Adopting A v = N H /{2.21 ± 0.09) x 
10 21 mag cm 2 jGiiver & Ozell l2009h and R v = 3.1, we find Ay = 0.4 8 ± 0.2 mag and 
E(B — V) = 0.15 ± 0.1 mag in agreement with lMason fc Cordoval (jl982al ). 



Figure [7] shows de-redden ed UV and V- and J -band fluxes calculated for E(B — V) = 
0.15 and the extinction law of ICardelli et al.l (119891 ) for Ry = 3.1. The top (green) line and 
triangles in the figure show the mean fluxes over orbital phases 0.75 - 0.85 and 0.15 - 0.25, 
during which the system is uneclipsed and contribution from the irradiated face of secondary 
star is minimized. The next (orange) line and squares show the fluxes during eclipse. The 
third (blue) line and crosses are the green minus orange fluxes, giving the flux from that 
part of the system that is eclipsed. All three of the SEDs roughly follow an F\ oc A~ 7 / 3 
distribution from the UV to the V-band, and then fall significantly below A" 7 / 3 in the J- 
band. The photometry of IMason fc Cordoval (Il982bl ) shows that the SED steepens further 
to a blackbody A~ 4 slope in the H and K bands. While it is tempting to identify the A -7 / 3 
SEDs with the theoretical SED for an optically-thick, steady-state a-model accretion disk, 
this identification is probably not correct. It is difficult to make an a-model disk produce 



2 Results provided by the ASM/RXTE teams at MIT and at the RXTE SOF and GOF at NASA's GSFC. 
http://xte.mit.edu 
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the same SED at mid-eclipse, when all but the outer parts of the disk are obscured, as it 
produces when it is not eclipsed and the entire disk is visible. Moreover, as we will show in 
the next section, most of the disk is hidden by optically-thick vertically-extended structures 
and is invisible at UV/optical wavelengths, so the theoretical SED of a simple a-model disk 
is not relevant. 

The bottom (red) line and open circles in Figure [7] show the mean fluxes at phases 
0.25 - 0.75 minus the uneclipsed (green) fluxes. This red SED isolates the flux in the large 
hump in the light curve at orbital phases opposite the eclipse. The flux contributed by 
the irradiated face of the secondary star is maximized at these phases. The flux from the 
irradiated secondary should be roughly a multi-temperature blackbody distribution. The red 
SED appears to be flattening near log A = 3.1, suggesting that the flux- weighted temperature 
of the heated face is ~ 25, 000 K. 

With the exception of the N V line, the flux in the UV emission lines is neither reduced 
during eclipse nor enhanced at other orbital phases. This shows clearly in the spectrum of 
the eclipsed parts of the system (the blue spectrum), which is nearly devoid of emission lines. 
The emission lines must, then, arise from optically-thin material far above and below the disk. 
After accounting for the extra width introduced by the poor spectral resolution of the UV 
spectrogram, we find the full width at zero intensity of the C IV emission line to be ~ 45 A. 
We compare this to its wid th in the HST Fain t Obj ect Spectrograph spectrograms from 
which the data presented by iPuchnarewicz et al. Jl995h were derived. In these spectrograms 



the width is near 40 A. These widths correspond to projected velocities up to ±4400 km s _1 . 
Thus the lines come from a high-velocity wind. Neither the strength nor the width of the 
C IV line depend on orbital phase. 

Since the material in the wind must be continuously replenished by outflowing gas closer 
to the disk, there must be a base to the wind. We do not observe C IV emission from the base 
of the wind, so either the base does not produce C IV emission or it is hidden by vertically- 
extended optically-thick material further out in the disk. The N V line is an exception. Its 
strength decreases by ~ 50% during eclipse and arises nearer the orbital plane. 
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5. Light Curve Analysis and a Model for 4U 1822-371 



5.1. Preliminaries 



5.1.1. The Light Curve Synthesis Program 

We modeled the UV/optical light curves of 4U 1822-371 usi ng a light curve synt hesis 
program^ we wrote to model X-ray binaries (e.g. UW Cor Bor, Mason et al. 2008 ) and 



related objects (e.g. SS Cygni, iBitner et al. 1 120071 ). The program assumes that the orbits 



of the stars are circular, the secondary star fills its Roche Lobe, and the primary star is a 
point source surrounded by an accretion disk. The accretion disk can have a complicated 
geometry. It can be non-circular, non-axisymmetric, and vertically extended; and it can have 
a rim, an interior torus, and bright or dark spots. The flux from the secondary star has a 
realistic spectrum using Kurucz stellar atmospheres (jKuruczlll996l ) but the remaining parts 
of the system emit like a blackbody. The program includes heating by irradiation. Because 
the geometry and temperature distributions in the model can be complex and asymmetric, 
the program resorts to ray tracing to calculate both irradiation and the orbital light curve. 
The program outputs synthetic light curves for Johnson/Cousins filter bandpasses and for 
square bandpasses over user-specified wavelength ranges. If provided with an observed light 
curve, our code calculates the x 2 f° r the fit of the synthetic light curve to the observed light 
curve. A calling program allows the user to find the model parameters that minimize \ 2 
either by a grid search through parameter space or by a simplex algorithm. 



5.1.2. The Masses and Dimensions of 4U 1822-371 

The arrival times of the X-ray pulses from the neutron star yield accurate values for 
the projected semi-major axis of the neutron star orbit a^s sin % = (3.015 ±0.015) x 10 10 cm 
and the mass function /(m 2 ) = m 2 g 2 (l + q)~ 2 siia 3 i = 0.0203 ± 0.0 003 where aMS is 



the radius of the neutron star's orbit and i is the orbital inclination (IJonker fc van der Klis 



200ll ). The orbital eccentricity is not measurably different from zero. 



It appears to us that the masses of the two stars in 4U 1822-371 have not yet been 
determined with certainty. The mass function can be recast to the form m^s sin 3 i = 
(0.0203 ± 0.0003)g~ 3 (l + q) 2 M . The eclipse solution limits the orbital inclination to a 
narrow range near 82°, so m^s ~ 0.0209g~ 3 (l + q) 2 M & . There have been several attempts 
to determine the mass ratio from the radial velocity curve of the secondary star. The spec- 



3 A full description of the program is available at http : //pisces . as . utexas . edu/robinson/XRbinary .pdf . 
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trum of the heated face of the secondary star has He I absorption lines and a narrow emission 



line f r om N III A4640, whi ch is excited by the Bowen fluorescence mechanism (IHarlaftis et al. 



19971 ; ICasares et al.ll2003l ). The radial velocity cu rve of the He I absorption has been mea- 
sured and its ampl itude is K HeI « 230 km s -1 (IHarlaftis et al.l Il997l ; ICowley et al.l 12003 



Casares et al.ll2003l ). It is, however, unclear how to interpret Ks e i because the absorption 
lines arise on only one face of the secondary and because the line profiles are probably dis- 
torted by helium emission lines from the accretion disk. The radial velocity curve also suffers 
from a phase shift of ~ 0.1P nrll with respect to the pulsa r ephemeris, further complicating 



its interpretation (IJonker et al.ll2003l ; ICasares et al.l 120031 ) 



The radial velocity curve of the N III A464 emission line also has been measured 

280 ± 3 km s" 1 ( Casares et al.l 120031 ; iMunoz-Darias et al- 



and has an amplitude K NT j 



20081 ). We agree with ICasares et al.l (120031 ) that this is a firm lower limit to the amplitude 
of the secondary's radial velocity curve, which translates to lower limits on the masses: 
nT'NS > 1-14 M & and m 2 > 0.36 M . Correcting from K^m to the true amplitude K 2 of 
secondary star's radial curve requires a model for the distribution of the N III A4640 emission 
across the face of the secondary star. IMunoz-Darias et al.l (120051 . 120081 ) have developed such 
a model and used it to deduce 1.52 M < m NS < 1.85 M & . In their model the N III A4640 
line is emitted uniformly over the face of the irradiated secondary except where shadowed 
by a simple thin disk with a raised, optically-thick rim. We will show that a more complex 
model is required for the disk in 4U 1822-371. The K-correction derived from the thin disk 
model is likely be only qualitatively correct. 

As the m asses of most neutron stars in binary X-ray pu lsars fall in a narrow range 
near 1.35M ( iThorsett fc Chakrabartyl Il999l ; iNice et al.l 120081 ). we adopted this value for 
most of our models. The mass function then becomes a relation between q and i, materially 
reducing the parameter space that needs to be explored when fitting models. This choice 
of neutron star mass also restricts other properties of 4U1822-371 to lie in a narrow range: 
If i lies in the range 78° < % < 84°, then 0.2948 < q < 0.3007, 0.398 < m 2 /M < 0.406, 
and the separation of the centers of mass of the two stars is 1.3336 < a/10 11 cm < 1.3367. 
Since the secondary star must fill its Roche lobe, its radiu s is close to 0.54 R^, pla cing it 
close to the observed main-sequence mass-radius relation ([Drilling fc Landoltl l2000l ). The 
temperature of a main-sequence star with this mass and radius is ~ 3500 K. This is a low 
enough temperature that we ignore flux from the secondary star except from those parts of 
its surface that are heated by irradiation. 



Because IMunoz-Darias et al.l (120081 ) argue for a high-mass neutron star in 4U 1822-371, 
we have also explored the effect of neutron star mass on our results, adopting the extreme 
value m^s = 2.0 M for a few models. For ttins = 2.0 M and i = 81° the mass function 
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yields q = 0.255, M2 = 0.51 M & and a = 1.50 x 10 cm. The radius of the secondary star is 
0.58 R Q , again placing it close to the main- sequence mass-radius relation. A main-sequence 
star with a mass of 0.51 M & has a temperature near ~ 3800 K. Fits of these higher-mass 
models to the light curves yielded small quantitative differences in the fitted parameters 
but the overall results were qualitatively almost identical. The derived orbital inclinations, 
for example, were typically higher by ~ 0.75° and the heights of the disk rim increased by 
10-15%. This also means, though, that the mass of the neutron star is not constrained by 
the optical/UV light curves. We report results only for the m NS = 1-35 M Q models. 



5.1.3. The Accretion Disk 

As discussed in the introduction, the unobscured X-ray luminosity of 4U 1822-371 
is high, possibly near Eddington; and the high mass flow rates deduced from the rapidly 
changing orbital period also suggest a luminosity near Eddington. Since the X-ray luminosity 
of 4U 1822-371 is relatively constant in the sense that no large-amplitude outbursts have 
been observed, the accretion disk in 4U 1822-371 is at least roughly in a steady state. These 
considerations suggest strongly that the accretion disk is optically thick and its inner edge 
is near the neutron star. 

While we will not necessarily limit our disk models to a-model disks, an optically-thick 
a-model disk is the default model for the outer parts of the disk in 4U 1822-371. The 
maximum radius of the disk is thought to be set by viscous dissipation induced by tidal 
interactions with the secondary star. For a mass ratio q = 0.3 the tidal truncation radius 



is r disk ~ 0.43a ~ 5.7 x 10 10 cm (IFrank et al.l 120021 ). Let us calculate a typical temperature 



and thickness of the disk at a radius that is large, but less than the radius of its outer 
edge. For a = 0.1, m^s = 1-35 M , a mass flow rate of 10~ 8 M yr" 1 , and a radius 
of r = 4 x 10 1 cm an optically- thick steady-state a-model disk has h/r^k = 0.04 and 



T e jf = 7800K ([Frank et al.l 120021 ) . These values depend only weakly on a, m ns , and m. 
Except where it is hit by the stream of transferred material from the secondary star, the 
outer, vertical rim of the disk should not be even as hot as this low temperature. We will 
see that the visible parts of the 4U 1822-371 system have a brightness temperature near 
3 x 10 4 K at UV/optical wavelengths. The temperature of the rim is so much lower that it 
is a minor contributor to the flux. If the rim is vertically extended, its inner surface may be 
heated to higher temperatures but this will only emphasize the low surface brightness of the 
outer surface of the rim. 



As a result, our models emit negligible flux from the disk rim. In the previously pub- 
lished models for the UV/optical light curve of 4U 1822-371, the disk rim is a significant 



13 



contributor, often th e dominant contr i butor to the UV/optical flux. This is shown most 
clearly in Figure 9 of iHellier &: Mason! (119891 ). Many of the differences between our results 
and earlier results flow from the different handling of the disk rim flux. 



5.2. The Disk Is Large, Its Center Obscured 

We first establish that the accretion disk is large and its center obscured at UV/optical 
wavelengths. The left panel in Figure [8] shows the eclipse portion of the binned UV light 
curve normalized to 1.0 just outside eclipse. The eclipse is wide, lasting A0 ps 0.20 in 
phase. The panel also shows two synthetic light curves, both calculated for a disk that is 
geometrically flat, has a nearly flat temperature distribution (T oc r -01 ), and is eclipsed at 
an orbital inclination % = 81°. One light curve was calculated for r^sk/a = 0.30, the other for 
?"disk/a = 0.45, and both have been normalized to match the flux at mid-eclipse to emphasize 
the effect of disk width on the light curve. Disks with r disk /a > 0.40 are needed to fit the 
observed eclipse. 

The dashed line in the right panel of Figure[8]is the synthetic light curve for a geometrically- 
thin, optically-thick, steady-state a-model disk with a T oc r~ 3//4 temperature distribution. 
The radius of the disk is r disk /a = 0.4, its temperature is 20,000 K at the outer edge, and the 
orbital inclination is 77°. The bright central regions of this model disk produce an eclipse 
light curve with a deep, narrow center in striking disagreement with the observed eclipse. 
The solid line in the panel is the synthetic light curve for a disk model with similar geometry 
but with a flatter temperature distribution to reduce the brightness of its central regions 
(T oc r~ ai and i = 78.5°). This model matches the general shape of the observed eclipse 
and, indeed, the fit is formally excellent. We do not conclude that the disk actually has a 
flat temperature distribution as we consider a flat distribution physically unrealistic. We do 
conclude that the flat temperature distribution is mimicking the true appearance of the disk, 
a disk that appears not to have bright central regions. As argued in the previous section, the 
central regions of the disk must be hot and bright, so the central regions of the disk must 
be obscured by disk structures extending above and below the orbital plane. 

We now consider two broad groups of models for obscuring the central disk. The first 
group includes models similar to standard thin-disk models. In these models the edge of the 
disk is tall enough to hide the center of the disk either because the disk has a large flare or 
because the disk is flat but has a vertically extended rim. In the second group of models the 
obscuration is produced by vertically-extended structures interior to the disk edge. 
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5.3. Obscuration by a Tall Disk Rim? 



The top panel of Figure [9] shows the observed UV eclipse light curve and a synthetic 
light curve produced by a geometrically-thin disk with a vertically-extended rim. The orbital 
inclination of the model is 81°, the disk radius is rdi S k/ a = 0.44, and the rim height is 
h r im/ a = 0.056, or h rim /r disk « 0.13. This rim is quite tall, raising concerns about the physical 
validity of the model; but the tall rim does produce a light curve that fits the observed light 
curve fairly well. Nevertheless, this model is fatally flawed. It has a combination of rim 
height and orbital inclination that leaves most of the central disk unobscured. Therefore, 
to avoid an eclipse that is too narrow and deep, the model must have a flat temperature 
distribution, thus failing to avoid the very problem it was contrived to solve. Worse, the 
neutron star is visible, disagreeing with X-ray observations. This is a generic and robust 
result for thin-disk tall-rim models. We have explored the parameter space for these models 
extensively, varying the orbital inclination, the disk radius and temperature, the rim height 
and temperature, and the irradiation of the inner surface of the rim by the inner disk and 
neutron star. Models that fit the observed light curve with an acceptable x 2 always have a 
flat temperature distribution and the neutron star is usually visible. 

Flared disks - disks whose thickness increase with radius - automatically have thick 
rims that can hide the inner disk and neutron star. In our computer code the height of a 
flared disk above the orbital plane is given by 



where a min and a max are the radii of the inner and outer edges of the disk, and .Hedge is 
the thickness of the disk at its outer edge. If H pow > 1, the disk has a concave flare. The 
middle panel of Figure M shows a synthetic light curve for a flared disk. The model disk has 
H e dge/ a = 0.060, a m i n /a = 0.01, a max /a = 0.43, and H pow = 2.0. The synthetic light curve 
fits the observed light curve about as well as the light curves for thin disks with thick rims 
but it has the same fatal flaws: The disk has a flat temperature distribution and the neutron 
star is usually visible. 

Models with concave flared disks have an important advantage over those with flat disks: 
The outer regions of the disk can be heated by irradiation from the neutron star and inner 
disk. The heating flattens the temperature distribution, providing a natural explanation for 
the flat distribution. The bottom panel of Figure M shows a synthetic eclipse light curve for 
a flared disk irradiated by the neutron star. The orbital inclination is 81°, a max /a = 0.45, 
Hedge/ a = 0.060, and H pow = 2.0. This value of H pow is much larger than the standard 
H P ow = 9/8 for an a-model disk, but lower values of H pow did not fit the observed light curve 
well. The fit for this particular model is good but the neutron star remains visible. 




(2) 
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In a final attempt to use the disk rim to block light from the center of the disk, we 
invoked a concave flared disk with an additional vertically-extended rim and allowed the 
disk and inner face of the rim to be irradiated by the neutron star. This model has enough 
free parameters to produce an excellent fit to the observed light curve at all orbital phases, 
not just the eclipse. But the inner disk and the neutron star still remain visible and the 
model fails. Thus, we reject this entire group of models. 

5.4. Torus Plus Variable-Height Rim 

The previous discussion leads us to consider vertically-extended disk structures closer 
to the neutron star. Inspired by the presence of the ADC at X-ray wavelengths, we give the 
disk vertical extension by adding an opaque torus (a donut with the neutron star centered 
in the hole of the donut). Our code allows the torus to have an elliptical cross section and 
specifies its geometry by three parameters: the distance from the center of the cross-section 
of the torus to the neutron star, which we call its radius or a a , the maximum height of the 
torus above and below the orbital plane, and the width of the torus between its inner and 
outer walls as measured in the orbital plane. The toroidal structure should not be taken too 
literally. Because 4U 1822-371 has a high orbital inclination, we observe the torus from the 
side and, like a donut viewed from the side, we observe only the outer wall of the torus. The 
observations give little information about the inner wall or about the interior of the torus. 
Furthermore, despite our original motivation for including the torus, we will see that the 
properties we derive for the torus are not consistent with the properties of an ADC and we 
will instead identify the torus with the optically-thick base of the disk wind. 

We assume that the torus emits blackbody radiation and for simplicity we give the 
surface of the torus a single temperature except where it is heated by irradiation. The 
temperature of the torus is poorly constrained. The SED of 4U 1822-371 is not a blackbody 
distribution, which precludes assignment of a meaningful color temperature. The large 
uncertainty on the distance and reddening to 4U 1822-371 preclude measurement of an 
accurate brightness temperature, but the brightness temperature of the outer wall of the 
torus must be surprisingly low. Fits to the UV eclipse show that the outer wall extends to 
~ 0.35a in the orbital plane (~ 80% of the disk radius) and the maximum height of the torus 
is ~ 0.2a (~ 50% of the disk radius), giving the torus a projected surface area of roughly 
0.28a 2 « 4 x 10 21 cm 2 . About half the projected area is eclipsed, A ed w 2 x 10 21 cm 2 , 
and the de-extincted eclipsed flux is f ec i = 3.1 x 10~ 14 erg/s/cm 2 /A at 1600 A (the blue 
line in Figure [7]). The brightness temperature of the eclipsed surface of the torus can be 
calculated from B x = f ed D 2 /A ech where D is the distance to 4U 1822-371. For D = 5000 pc 
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the monochromatic brightness temperature at 1600 A is T5 m 26, 000 K. Although this 
temperature cannot be trusted to within a factor of two or so, it is ~ 10 3 times lower than 
temperatures typically associated with coronae. Thus, the brightness temperature of the 
torus - at least the outer wall of the torus - is not consistent with an ADC. The temperature 
is, however, high enough that the torus dominates the UV/optical flux from 4U 1822-371. 
This justifies our earlier claim that the flux from unirradiated surfaces of the secondary star 
and the outer surface of the disk rim can be ignored. 

This torus-like structure is far too tall to agree with the thickness of an a-model disk 
for any reasonable choice of parameters. Its temperature is far too low for it to be part of 
the ADC. What, then, is this vertical structure? Since we have already found a high velocity 
disk wind and have shown that there is much mass loss from the disk, we suggest that the 
vertical structure is simply the optically-thick base of the disk wind. 

The torus by itself does not produce a synthetic light curve that fits the observed eclipse 
light curve adequately, but the addition of either an opaque vertically-extended rim or a flared 
disk yields an excellent fit. The best-fit synthetic light curve for a flared disk is shown in 
Figure [TU1 The model from which the light curve was calculated has an orbital inclination 
of 83.5°, a concave flared disk with radius r disk/a = 0.43 and edge height -f/edge/^disk = 0.03. 
The midpoint of the torus is at a Q = 0.53rdi S k, its maximum height is 0.56rdi S k, and its width 
is 0.56rdisk, so that the torus extends from 0.26rdi S k to 0.81rdi S k- The torus blocks flux from 
the neutron star and all the disk out to 0.81rdi S k, and the edge of the flared disk blocks 
flux from most of the rest of the disk. Thus the torus solves the two problems that vitiated 
models with flat or flared disks and raised rims. 



5.5. A Model for the Entire Orbital Light Curve 

We now model the entire orbital light curve using the eclipse model from Section 5.4 
as a starting point. The salient feature of the non-eclipse part of the orbital light curve is 
the large-amplitude, roughly-sinusoidal hump. The phase of the hump might lead one to 
suspect that the hump is caused primarily by the varying visibility of the irradiated face of 
the secondary star but our attempts to model the hump with just the irradiated secondary 
star yielded poor fits. The hump's precise phase, its shape, and the variations of its shape 
from year to year cannot be produced by the irradiated secondary alone. 

The model producing the synthetic light curve shown in Figure (TU] had a torus and a 
disk with a dark, vertically-extended edge. We now allow the height of the edge to vary with 
azimuth around the neutron star. As its height varies, the rim blocks more or less of the 
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flux from the torus, helping to produce the hump. In the previous section we introduced the 
vertically-extended edge by invoking a concave flared disk, but a flat disk with a vertically 
extended rim would have worked just as well. In the current context it is more convenient 
to model the edge with a flat disk and extended rim, allowing the rim height to vary with 
azimuth. 

We also include flux from the heated secondary in the model and allow fits to the light 
curve to determine how much flux the secondary contributes. The secondary is heated in part 
by the torus and disk but mostly by X-ray flux. While we include X-ray flux in the model, it 
is not a physically meaningful quantity since the amount of flux needed depends sensitively 
on the geometry of the X-ray emitting regions, on the radiative transfer of the flux through 
the ADC and torus, and on the reprocessing and redistribution of the flux in the secondary 
star. As we have little information about any of these factors, we regard the X-ray flux as 
merely a convenient parameter to specify the irradiative heating of the secondary star. With 
the addition of a variable-height rim and flux from the irradiatively-heated secondary star 
we achieved successful fits to the entire orbital light curve and no further complications were 
needed. 

While the foregoing discussion has been necessarily lengthy, the resulting model is ac- 
tually rather simple. The parameters describing the geometry of the model are 

• The orbital inclination, which also fixes the mass ratio and the mass of the secondary 
star through the mass function, 

• the outer radius of the disk, 

• the height, width, and radius of the disk torus, and 

• the height of the disk rim as a function of azimuth around the neutron star. 

We assume that the secondary star and the disk rim have low temperatures and do not 
contribute significant flux except where heated by irradiation. Fits to the light curve show 
that the inner regions of the disk are hidden by the torus and the outer regions are hidden 
by the disk rim, so the disk is not visible and its temperature is irrelevant. Thus the only 
temperature that needs to be specified is the monochromatic brightness temperature of the 
torus. Finally, the amount of X-ray energy irradiating the disk rim and the secondary star 
must be specified. There are also a few nuisance parameters such as the zero point in orbital 
phase, a scale factor on the flux, and various albedos that must be specified but are irrelevant 
to the model. 

We fit the model to the observed UV light curve using a simplex algorithm to find 
the set of parameters that minimized the \ 2 °f the fit. We also examined parameter space 
extensively with grid searches to test the uniqueness of the fit. Within the context of the 
model the values of the fitted parameters are unique and robust. The resulting synthetic 
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light curve is superimposed on the observed light curve in the top panel of Figure [HJ The 
rim height as a function of azimuth angle is shown in Figure [121 where the azimuth angle 
is defined to be 0° on the side of the disk opposite the secondary star and increases in the 
direction opposite the direction of orbital motion. We fit only the V- and J-band light curves 
obtained in 2006, the same year as the UV light curves were obtained. To fit the light curves 
we retained the values of all the parameters determined from the fit to the UV light curve 
except for those specifying the geometry and temperature of the torus. The fitted synthetic 
light curves in the V- and J-bands are shown superimposed on the observed light curves in 
the bottom two panels of Figure[TlJ The values of x 2 f° r the fits are 163.5 for 195 data points 
and thirteen free parameters in UV, 210.3 for 150 data points and four free parameters in V, 
and 130.7 for 116 data points and four free parameters in J, for a total reduced \ 2 of 1.05. 

Values for the fitted parameters are given in Table [2J The radius of the disk is r^k/a = 
0.40 ± 0.01 and the disk rim can have a height up to /i r i m /^disk ~ 0.17, but the height for 
much of the disk is near /WiA'disk ~ 0.05, close to the a-model h/r values (Section 5.1). The 
height of the torus, 0.20a = 0.5rdi S k, is the same in all three passbands, but the position of 
its outer wall is much different, a Q = 0.31a = 0.78rdi S k in the UV, but a Q = 0.39a = 0.98ra s k 
in the V- and J-bands. This is a reflection of the obvious difference between the UV eclipse 
and the broader, flatter bottomed eclipse at longer wavelengths. 

This model also fits the 2005 V- and J-band light curves, but with a change to the 
rim height with disk azimuth. The maximum disk rim height is approximately the same 
as the 2006 model with h rim /r disk ~ 0.15. However, the rim is now symmetric about the 
approximate location of the "hot spot" where the tidal stream encounters the disk and not 
the radius connecting the companion star and the neutron star. A detailed quantitative 
model is unwarranted due to the sparse sampling of the data. 



On considering these results, a new issue arises. While we and iHellier fc Mason! (119891 ) 
before us needed to place the opaque wall near the edge of the accretion disk to fit the light 
curve, neither the X-ray nor the UV/optical light curves require it to be precisely at the disk 
edge. It could, in fact, be anywhere between the edge of the disk at ~ 0.4a and the outer 
wall of the torus in the V- and J-bands at ~ 0.39a. Although this is not a large change in 
position, it allows a significant re-interpretation of the wall. We placed the wall at the disk 
rim - and we suspect others placed it at the disk rim - primarily because there was no other 
disk structure to which the wall could be attached. With the addition of a disk torus the wall 
can plausibly be made part of the torus, becoming an opaque belt around the outside of the 
torus. The wall now becomes part of the disk wind, not a disk rim. This re-interpretation 
has a considerable advantage. When made part of the disk wind, the wall no longer need be 
in hydrostatic equilibrium, requiring impossibly high gas pressure or supersonic turbulence 
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for support. It is simply a cooler, darker layer at the base of the disk wind. 

6. Summary and Discussion 

Combining our new times of eclipses with the previously published times, we have 
derived an improved ephemeris for the eclipses in the UV/optical light curves of 4U 1822- 
371 (equation [Q . The quadratic term in the ephemeris yields a timescale for a change in 
orbital period of P/P = (3.0 ± 0.3) x 10 6 yr. The deduced rate of mass transfer is large, 
probably greater than 3 x 10~ 8 M yr -1 . Mass cannot be accreting onto the neutron star 
at this rate without violating the Eddington limit, so much of the transferred mass must be 
lost from the system. 

We have inferred a strong disk wind from the C IV emission line in the UV spectrum. 
The wind is axially symmetric and has projected outflow velocities up to 4000 km s^ 1 . The 
flux in the C IV emission lines does not decrease during eclipses. Since regions less than 
~ 0.5rdisk above the disk are at least partially eclipsed, the C IV emission arises from regions 
yet further above the disk. The outflowing material closer to the disk either does not produce 
C IV emission or is hidden by other optically-thick structures. 

Much of this paper was devoted to developing a model for the accretion disk in 4U 1822- 
371 from fits to the UV/optical orbital light curves. The following properties of the accretion 
disk appear to be robustly determined. The disk is large; the fitted radius is r^k/a = 
0.40 ± 0.01, close to the tidal truncation radius. To avoid an eclipse that is deeper and 
narrower than the observed eclipse, the neutron star and the central regions of the accretion 
disk must be obscured. The dual requirements that the models must produce synthetic light 
curves that fit the observed eclipse and must obscure the neutron star and central disk forces 
most of the obscuring material to be at intermediate radii in the disk, not at the edge of 
the disk. We modeled the obscuring material as an optically thick torus. The torus extends 
more than 3/4 of the way to the edge of the disk at UV wavelengths and nearly all the way 
to the edge in the V- and J-bands; and it extends to a height of ~ 0.5r disk above the disk at 
all wavelengths. 

The surface brightness of the outer (visible) wall of the torus is only ~ 26, 000 K. 
While this estimate is highly uncertain, the true brightness temperature is unlikely to be 
more than a factor of 2 - 3 higher and is far less than coronal temperatures. We have, 
therefore, interpreted the torus as the optically-thick base of the disk wind, not as the 
optical counterpart of the X-ray ADC. 

We modeled the large hump in the orbital light curve of 4U 1822-371 with a combination 
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of heated face of the secondary star and a variable- height, optically-thick wall located at the 
edge of the disk. All previous models for the orbital light curve have needed a similar wall, 
but we now suggest that the wall is actually slightly closer to the neutron star than the edge 
of the disk and is part of the disk wind, not the disk rim. The wall is a cooler, darker layer 
at the base of the disk wind just above the disk. 

Our results thus yield a more complex view of the structure of the disk in 4U 1822-371. 
The previous X-ray observations have shown that there is a vertically-extended, optically- 
thick ADC stretching roughly half way out to the edge of the disk. Our UV/optical obser- 
vations show that most of the outer disk emits a cooler disk wind. The wind is is optically 
thick close to the accretion disk but becomes optically thin at a height of ~ 0.5rdi S k above 
the disk. Close to the disk the wind is yet cooler and darker, forming a relatively dark wall 
around disk. The disk itself appears to be nearly entirely buried within and obscured by 
these vertically-extended structures. 

We thank Raymundo Baptista and Coel Hellier for their valuable insights into this work. 
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Fig. 1. — The UV spectrum of 4U 1822-371. The prominent emission lines are N V at 
1240 A, a blend of O IV, O V, and Si IV near 1370 A, the C IV doublet at 1548 A/1550 A, 
and He II at 1640 A. There are also ISM absorption blends of Si II/C I and O I/Si III near 
1260 A and 1300 A respectively. There is a significant red leak in the SBC at wavelengths 
longer than ~ 1900 A. 
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Fig. 2. — The O — C diagram for the optical/UV eclipse ephemeris (equation [[J. The open 
circles are the V eclipse times from SMARTS data and the open squares are the UV eclipse 
times from HST data. The solid circles are the previously published times of optical eclipses 
listed in Table [TJ The arrows show the mean offset between the ephemerides for the X-ray 
and UV/optical eclipses. 
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Fig. 3. — The UV light curve of 4U 1822-371 from 2006. The fluxes were extracted from the 
HST spectrograms and folded on the orbital period of 4U 1822-371 (equation [TJ to produce 
one light curve for each of the HST visits. The internal errors on the fluxes are ~ 0.5%, 
giving error bars smaller than the points used to plot the data. 
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Fig. A.— V-band light curves of 4U 1822-371 in 2005 and 2006. The data have been 
folded on the orbital period given by the UV/optical ephemeris (equation [T]). The error bars 
correspond to the internal measurement errors. 



-27- 




Fig. 5.— J-band light curves of 4U 1822-371 in 2005 and 2006. The data have been 
folded on the orbital period given by the UV/optical ephemeris (equation [1]). The error bars 
correspond to the internal measurement errors. 
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Fig. 6.— The X-ray light curve of 4U 1822-371 between 1996 January 5 and 2008 March 
5 from the All Sky Monitor on RXTE. Each point is the monthly average of all three ASM 
bandpasses. The vertical dashes mark the times of e clipse minima at U V/optical wave l ength s 
measured by I Cowley et al.l (120031 ) , iBaptista et al.l (120021 ) (4 eclipses) , ICasares et al.l (120031 ) , 
and by us with EST and SMARTS. 
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Fig. 7. — The dereddened UV and V- and J-band fluxes. The top (green) line and triangles 
show the mean fluxes at orbital phases 0.75 - 0.85 and 0.15 - 0.25, during which the system 
is uneclipsed and contribution from the irradiated face of the secondary star is minimized. 
The next (orange) line and squares show the fluxes during eclipse. The third (blue) line 
and crosses are the green minus orange fluxes, giving the flux from that part of the system 
that is eclipsed. The bottom (red) line and open circles show the mean fluxes at phases 
0.25 - 0.75 minus the uneclipsed (green) fluxes, which should be dominated by the flux from 
the irradiated face of the secondary. The solid line shows fx oc A -7 / 3 and the dashed line 
fx oc A" 4 . 
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Fig. 8. — The solid circles show the co-added and binned UV eclipse light curve of 4U 1822- 
371 normalized to 1.0 just outside the eclipse. Left Panel: The lines are synthetic light curves 
for a geometrically-flat disk with a flat temperature distribution (T oc r -01 ) and for an orbital 
inclination of 81°. The solid and dashed lines are the light curves for disks with r disk/a = 0.45 
and 0.30 respectively; both have been normalized to emphasize the eclipse shape. The disk 
in 4U 1822-371 must be large, rdi S k/ a ~ 0.40, with its precise radius depending on the orbital 
inclination and other disk properties. Right Panel: The dashed line is the eclipse light curve 
for a geometrically-flat disk with an a-model temperature distribution (T oc r~ 3 / 4 ) and a 
surface temperature of 20,000 K at the maximum radius, but 2,000 K at the outer edge. 
The solid line is the synthetic light curve for a disk with the same surface temperature at 
its maximum radius but a nearly flat temperature distribution (T oc r~ ai ). The synthetic 
light curves have been normalized to emphasize the shape of the eclipses. The a-model disk 
produces an eclipse that is much too narrow to fit the observed eclipse. 
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Fig. 9. — The solid circles in all three panels show the co-added and binned UV eclipse 
light curve of 4U 1822-371 normalized to 1.0 just outside the eclipse. The solid line in 
each panel is a synthetic light curve fitted to the observed eclipse. Although the synthetic 
light curves fit the data well, they are all fatally flawed. The central regions of the disk are 
unobscured, so the models must have a flat temperature distribution, thus failing to avoid 
the very problem they were contrived to solve. The neutron star is also visible, disagreeing 
with X-ray observations. Top Panel: The synthetic light curve for a thin-disk tall-rim 
model. The orbital inclination of the model is 81°, the disk radius is r^^/a = 0.44, and the 
rim height height is hn m /a = 0.056. Middle Panel: The synthetic light curve for a concave 
flared disk. The orbital inclination of the model is 78°, a max /a = 0.43, H edge /a = 0.060, and 
H P ow = 2.0 (see equation [2]). Bottom Panel: The synthetic light curve for a concave, flared, 
irradiated disk. The orbital inclination of the model is 81°, a max /a = 0.45, H e ^ e /a = 0.060, 
and H pow = 2.0. The disk is irradiated by flux from the neutron star. 



-32 - 




o 



-0.1 



0.1 



Orbital Phase 

Fig. 10. — The solid circles show the observed UV eclipse light curve of 4U 1822-371. The 
solid line is a synthetic light curve fitted to the data. The model from which the light 
curve was calculated has an orbital inclination of 83.5°, a concave flared disk with radius 
r disk/ a = 0.43 and an edge height -ffedgcAdisk = 0.07. The disk is covered by an opaque torus 
with midpoint at a Q = 0.53rdi S k, maximum height of 0.56rdi S k, and width of 0.56rdi S k, so that 
the inner wall of the torus extends in to 0.26rdi S k and the outer wall out to 0.81rdi S k- The 
torus blocks flux from the neutron star and from the disk out to 0.81r d i S k! and the edge of 
the flared disk blocks flux from most of the rest of the disk. 
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Fig. 11. — The best fitting UV, V, and J light curves with the best fitting model in Table 
[2J This model includes an optically thick torus. 
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Fig. 12. — The variation in disk rim height with azimuthal disk angle. An angle of 0° is away 
from the star, rotating opposite the direction of orbital motion, making 180° the radius in 
the disk connecting the NS and secondary star. The rim height is in terms of disk radius 
and is the height added to the very thin (0.002r disk ) disk. 
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Table 1: Times of T min and 0-C values as calculated using a 
quadratic least squares fit to the ephemeris. 



Eclipse Time (HJD) 


Error in 


Ref. 


0-C 


2444044 8450 


0060 

\J • \J \J V7 V7 


1 


-0.0207 


2444090.1140 


0080 

V 7 > V 7 V 7 v 7 V 7 


1 


-0.0128 


2444101.0280 


0.0080 


1 


-0.0079 


2444105 6650 


0060 

V7 • V7 V7 V7 V7 


1 


-0.0131 


2444106 5970 


0060 

V7 • V7 V7 V7 V7 


1 


-0 0095 

V7 • V7 V7 t7 


2444137 9350 


0.0100 


1 


-0 0062 

V7 • V7 V7 \J 


2444411.1320 


0080 


2 


-0.0007 


2444412 0580 


0080 

V7 • V7 V7 vJ V7 


2 


-0.0031 


2444664 8120 

J_ _l_ _l_ V./ V7 -L • ^7 _L V 7 


0080 

V7 • V7 V7 V7 V7 


3 


-0.0152 


2444783 8940 


0060 

V7 • V7 V7 V7 V7 


2 


-0.0048 


2445579 5650 

J- J- ..J .J 1 » ./ • / V7 .J V 7 


0030 


4 


-0.0019 


2445580 7250 


0030 

V7 • V7 V7 KJ V7 


4 


-0 0025 

V7 • V7 V7 i—1 


2445615 3115 

£*J _L _L v../ V./ -L ".7 > X _L '7 


0020 

V7 • V7 V7 V7 


4 


-0 0002 

V7 • V7 V7 V7 ^ 


2445937 7110 

fcJ Al Al '.J kJ -.7 1 • 1 J L \J 


0.0021 


4 


0.0004 


2446234 5787 


0.0018 


4 


0.0010 


2447296 4798 


0.0010 


4 


0.0040 


2447379 3410 

JI JI I (J 1 ^-7 JI X_ V/ 


0.0010 


4 


0022 


2447999 7674 

^— JI JI | (7i7 ■ 1 V./ 1 IT 


0039 


5 


0.0011 


2449163 0960 


0028 

V7 • V7 V7 i-J vJ 


6 


-0 0020 


2449164 0237 


0028 

V7 • V7 V7 Vj 


6 


-0.0027 


2449164.2542 


0028 

V 7 ■ V 7 V 7 Y7 


6 


-0.0043 


2449165.1852 


0028 

V 7 • V 7 V 7 V7 


6 


-0.0018 


2449166 1190 

i—i _l_ _l_ t.7 _L V7 V7 • 1 1 r.7 V 7 


0028 


6 


0036 

V7 • V7 V7 \J 


2450250 5308 


0008 

V7 • V7 V7 V7 W 


7 


-0 0003 

V7 • V7 V7 V7 


2450250 7626 

i—i _l_ v.7 V 7 v.7 V 7 > 1 V7 V 7 


0008 

V7 • V7 V7 V7 W 


7 


-0 0006 

V7 • V7 V7 V7 V7 




n nnns 

U.UUUo 


7 


n nnn/i 

-U.UUU4 


2451264.8446 


0.0035 


5 


-0.0060 


2451373.7094 


0.0008 


7 


-0.0007 


2451756.4577 


0.0012 


7 


-0.0016 


2451782.4542 


0.0012 


7 


-0.0015 


2452089.7692 


0.0008 


7 


-0.0001 


2453618.6767 


0.0050 


8 


-0.0023 


2453828.9720 


0.0010 


8 


0.0010 


2453830.8299 


0.0010 


8 


0.0020 


2453932.7201 


0.0040 


8 


-0.0042 



(1) 

(5) 



Mason et alj Jl980l) (2) iMason et al.1 (1 19821). ( 3) ICqwlev et all ll 19821). tt) iHellier k Maso3 (| 19891) 



Cowlev et all |2003l) . (6) lHarlaftis et al.l (jl997l ). (7) lBaptista et all j20od) . (8) This paper 
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Table 2: Parameters of 4U 1822-371 as determined by model light curves. 



Parameter 


Value 


Period (days) 


0.232109 


Inclination (degrees) 


81.17 ± 0.14 


M NS (M ) (fixed) 


1.35 


q = M 2 /M NS 


0.2971 


-^pow 


1.1 


' aisK 


0.4 ±0.01 


UV Torus 




Location, ao 1 


0.165 


Torus Width 1 


0.283 


Torus Edge Height 1 


0.201 


V Torus 




Location, ao 1 


0.245 


Torus Width 1 


0.282 


Torus Edge Height 1 


0.205 


J Torus 




Location, ao 1 


0.243 


Torus Width 1 


0.288 


Torus Edge Height 1 


0.204 



x In units of the orbital separation (see Section 5.1.2). 



